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PACS 75 . 45 . + j - Macroscopic quantum phenomena in magnetic systems 

Abstract. - We study the quantum dynamics of polygonal rare-earth molecular clusters with 
Ising-type ion magnetization. It is shown that the ground state of such systems is a non-magnetic 
quasi-doublet of states with oppositely twisted ion spins. The states differ in sign of toroidal 
moment, which is a natural physical quantity to characterize the spin chirality of the clusters. The 
possibility of macroscopic quantum tunneling of toroidal moment between the states is predicted. 
The effects of an external current is considered, both in equilibrium and in the frames of the 
Landau-Zener-Stiickelberg tunneling model. The special treatment is given for the most important 
case of triangular rare-earth clusters. 



Introduction. — Magnetic nanoclusters are in the fo- 
cus of intense research because they exhibit remarkable 
quantum properties such as macroscopic quantum tunnel- 
ing (MQT) of magnetization, molecular bistability, quan- 
tum interference, Berry phase effects, etc. Besides the ob- 
vious fundamental significance, their investigation is of 
great importance for the fronts of nanoclcctronics, namely 
molecular spintronics [l]-[4] , as well as for full-scale quan- 
tum computer development in the aspect of overcoming of 
non-scalability and difficulties in qubit state control. 

In this respect, magnetic molecular nanoclusters with 
spin chirality (Cu 3 , V15, Dy 3 , etc.) are of special interest. 
The extra degree of freedom can be used for qubit cod- 
ing. It was shown by the example of Cu 3 magnetic cluster 
[5] that the chirality is the source of a spin-electric effect 
whereby the spin states of the cluster is driven by an elec- 
tric field, which is better localized than magnetic one |6j. 

The alternative possibility is the clusters with the non- 
magnetic ground state. It is pointed [7| that getting rid of 
the spin makes the qubit very little sensitive to magnetic 
noise, but creates at the same time a potential problem 
for measurement since there is no Zeeman coupling of the 
qubit to an external magnetic field. One of possible solu- 
tions to this problem relies on the presence of an orbital 
moment associated to the chirality [7]. 

Recent experimental and theoretical studies [8 -10 of a 
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dysprosium based triangular cluster demonstrate the pos- 
sibility of chirality control with an electric current (or just 
with crossed electric and magnetic fields) via interaction 
with toroidal moment, which is a natural characteristic of 
chirality 
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It is remarkable that the cluster is of zero 
magnetic moment in the ground state in spite of the fact 
that there is a certain order in spin arrangement in it. The 
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molecules with toroidal moment were first implied in 
Recently, such a molecule (namely Dy 3 ) was synthesized 
12 and experimentally investigated [8](9]. Other examples 



of molecules with toroidal moment are nanocluster V15 
13 and probably nanocluster Cu 3 j5j. Toroidal ordering 



in crystals is reviewed in 14 



Today, the quantum properties of chiral molecular clus- 
ters with toroidal moment is lacking of systematic re- 
search. In the present paper, the intriguing question of 
macroscopic quantum tunneling of toroidal moment (spin 
chirality) is considered for the first time. We predict the 
existence of Rabi oscillations in the clusters, which pro- 
vides the good evidence for the MQT. The current driven 
dynamics is also considered, both in equilibrium and in the 
frames of the Landau-Zener-Stiickelberg tunneling model. 
Due to the coupling between the toroidal moment and the 
current, it is possible to distinguish between states with 
opposite chirality, which is important for the purpose of 
qubit coding. 

The model. — We consider a spin ring, a system of N 
non-Kramers rare-earth ions located in the apices of a reg- 
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Fig. 1: The spin structure of a triangular rare-earth molecular 
cluster and the local easy axes orientation in respect of the 
laboratory XYZ reference frame. The thicker arrows represent 
the spins of the rare-earth ions in the molecule in the state 
with toroidal moment Tz = +3Tq. 



ular polygon. The Hamiltonian of the rare-earth polygonal 
cluster reads as follows 



H 



N 



N 



Hint + 9jI^b ^ HJ,, 



(1) 
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where is the operator of the crystal field at the i-th 

ion location, is the total angular momentum of the z-th 
ion, H is the external magnetic field strength, and Hint 
is the Hamiltonian of the dipole and exchange interactions 
of the rare-earth ions in the cluster. 

The ground state of rare-earth ions in compounds is 
formed mainly due to the influence of a crystal field and 
often has a strong anisotropy of the magnetic moment. For 
example, the ground state of Dy 3+ ions in Dy 3 cluster is 
very close Kramers doublet \Mj — ±15/2) [8}{l0] and re- 
sponds only to the z% local component of an external mag- 
netic field (see fig. 1). The first excited state is separated 
from the ground one by the energy of 200 cm -1 . The wave 
functions of the excited state are close to \Mj — ±13/2). 
We can therefore conclude that the environment of rare- 
earth ions is almost axially symmetric. In the case of non- 
Kramers ions, small asymmetrical perturbations remove 
degeneration, thus making the ground state be quasi- 
doublet, i.e. the two close singlets with splitting small com- 
pared with distance from the higher levels. According to 
the Griffiths theorem 
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the magnetic moments of the 
ions can be directed only along a specific axis, namely the 
local Zi-axes. 

Let |cij) and |6j) be the eigenstates of the z-th (hereafter 
i = 1, N) ion in the crystal field in respect of the local 
axes. The z^-axes are perpendicular to the bisectors of the 
iV-gon. The abscissa then directed along the 

bisectors at the apices of the A-gon. Projection of the 
Hamiltonian in eq. ([TJ on the Hilbert space with functions 

\x± ) = (\o-i) ± i\bi))/V2 constituting the basic set yields 
the effective Hamiltonian 
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where a x , a y , and a z stand for the Pauli matrices, j is 
the exchange interaction constant, A is the energy gap 
between the singlet levels \ai) and \bi) (as a rule A << 
j), and fii = gjn B Im(ai\J Zi \bi). Usually, \x ( ±) ~ \Mj = 
±J) = |±). It will be shown below that eq. ^ is valid not 
only for non-Kramers ions, but for Kramers ions as well. 
In the latter case, the splitting A can be produced by an 
external magnetic field H, and A — > at H — >• 0. 

There are 2 N possible orderings of the ion spins in the 
rare-earth polygonal cluster. The wave functions of the 
cluster can be written as 
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where a ni stands for the sign ("+" or "— ") of the i-th ion 
(i = 1, A) spin projection onto the local Zi axis in the 
ra-th state (n = l,...,2 N ). 

The spin orderings can be characterized in terms of 
spin chirality. It is clear that spins in the states \n (n = 



1,...,2 N and their conjugates \ n (n = 2 
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are inversely twisted, i.e. the states have opposite chirality. 
The natural physical quantity associated with spin chiral- 
ity in this case is the T-odd polar vector of the anapole 
(toroidal) moment, which corresponds to the first term of 
the toroidal family in the multipole expansion of an arbi- 
trary electric current distribution. We would remind that 
the toroidal moment operator in the case of localized mag- 
netic ions can be defined as 



N 

9jVB^2[ri X Ji], 

1=1 



where is the radius-vector connecting the center of the 
A-gon with the i-th apex (all r, = ro). 

Hereafter we deal with the dimensionless values of the 
toroidal moment r = Tz/T , where T = ^jUb^oJ ■ Ob- 
viously, the toroidal moment of the system in a Xn state 
is t = q + — <7_, where q± is the number of "+"s and "— "s 
in the Xn state. The values of the toroidal moments for 
conjugate states are different in sign. 

Macroscopic quantum tunneling of the toroidal 
moment. — In the limiting case A —> (i.e. A << 
j) the ground state is | + ...+) and | — ...—), which is 
doubly degenerated. Such a situation brings up the actual 
question on the macroscopic quantum coherence (MQC), 
i.e. the cluster ground state degeneration removal, and on 
the macroscopic quantum tunneling (MQT) of the toroidal 
moment, for instance, from the state \Tz = +T) to the 
state \T Z = -T), where T = NT a + 0{(A/j) 2 ) is the 
value of the toroidal moment at the finite quantity A. Both 
questions come to calculation of the imaginary-time (r = 
it, t £ [0;to]) transition amplitude between the states, 
which can be treated in terms of the path integral |16| 
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where DODip is the integrating measure, 



/•To 

Se = C{— ir)dr 

Jo 



is the Euclidean action, and C is the Lagrangian of the 
system. The wave functions of the rare-earth polygonal 
cluster in the considered Hilbert space, which is the direct 
product of the state spaces of each of the N cluster ions, 
can be presented as 



N 



(4) 



where a(9,tp) = cos | and /3(9,tp) = sin | • e~' lip are the 
Cayley-Klein parameters depending on polar 9i and az- 
imuthal tfi angles of the <7j in the i-th local coordinate 
system with polar axis directed along the i-th easy axis. 

The quantum dynamical equation of motion for the 
toroidal moment reads 



ih- 



dT 



[T,H], 



(5) 



where the Hamilton operator % is defined by eq. We 
suppose that the dynamics of T(t) is mainly conditional 
on the spins of the cluster ions, i.e. rj(i) = const. If one 
averages out eq. ([5j with the wave function from eq. Q, 
one obtains 



2~<tt 



rii x 



dE 
dm 



(6) 



where = {sin#i cos tpf, sinf^ sin tpf, cos are the unit 
vectors and 



A 



N 



(7) 



is the energy of the system at the zero temperature. 

It is crucial for the MQT to take place that there are co- 
herent macroscopic processes, which naturally means that 
9x = ... = 6n = 9 and tpi = ... = tp N = tp. The symmetry 
of the system is not affected unless the external magnetic 
field has a component directed in the plane of the rare- 
earth polygon. In terms of the (6, tp) variable set eqs. ([6| 
take the form 



h(p = 4j cos 9 — A cot 9 cos tp, 
h9 = — A sin<^>, 



(8) 



which are the Lagrange-Euler equations of the Lagrangian 
C = ^(1 - cos9)<p + Nj cos 2 9 + sin 9 cos tp. (9) 

The first term in eq. ^ is the Wess-Zumino one. Gen- 
erally, the term is calculated making use of well-known 
formula for the Berry phase and interpreted as a contribu- 
tion conditional on the gauge field of the Dirac monopole 
located at the point n = (0; 0; 0). It is interesting that the 



Isingness of the considered system do not have an impact 
on the form of the Wess-Zumino term, i.e. on the interac- 
tion between the molecule and the geometry field of the 
Dirac monopole, as opposed to its interaction with the 
external field. 

The probability of the tunneling per time unit in the 
quasi-classical limit is determined |17| by the expression 

P = Aexp(-B), 

where B = 2Se/Ti is the Gamov constant, and 

S E = J (-i^(l-cos9)^- + E(9,tp))dT 

is the Euclidean action calculated on the instanton tun- 
neling trajectory. The pre-exponential factor A ~ ujq 
where ujq is the instanton frequency (see below). 

In order to calculate the tunneling trajectory and the 
corresponding contribution to the action we switch ana- 
lytically in eq. (|8| and eq. ^ to the imaginary time r = it, 
and note that eqs. ^ have the first integral 
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sin 9 cos tp = const. 



(10) 



Now we consider the stationary points and the transi- 
tion solution ("bounce trajectory") of the eqs. Points 
tp = and cos 9 = ±00 = iyl — (A/4j) 2 are the station- 
ary points of the system, which are actually the terminal 
points of the tunneling between states with toroidal mo- 
ment T z =±T = ±NT ■ a . Point 9 = ir/2 and tp = is 
the minimum point of the inverted energy. 

Eliminating the ^-variable from the system of eqs. ([8]), 
we come to the instanton equation 



dcos9 
dr 



2 1 
Ti 



(11) 



Integrating the equation and taking into account the ini- 
tial conditions cos0| r= ±oo = ±00, we get the instanton 
solution 

a , ir, f ^ot" 

cos 9 = do • tanh — — 

V h 

whence it follows that the instanton frequency lu is 
The time dependence of the tp- variable is then specified by 



the motion integral, see eq. ( 10 ) 
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It should be noted that Euclidean angle tp takes on imagi- 
nary values. Substituting the obtained explicit dependen- 
cies 9 = 9(t) and tp — tp{r) into eq. ^ for the system 
Lagrangian and calculating the action Se on the tunnel- 
ing trajectory, we obtain for the Gamov constant 
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dr, (12) 
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which finally comes to 



B = N In 



1 + O-Q 
1 - (T 



2a 



(13) 



It is seen that B = 2A<Tq/3 if Oq — > 0, and also B — >• oo 



if <To — > 1. The tunneling frequency is f 



\ It is 



determined by the <Jo, the value of the effective toroidal 
moment. If the exchange interaction between the cluster 
ions is strong enough (j >> A), then <tq — >• 1, and the 
tunneling probability vanishes. 

The quantum mechanism of the tunneling dominates 
over thermal-activated processes of the spin reorientation 
at low enough temperatures. The temperature dependence 
of the speed of the latter is ~ exp(— Uo/kT), where Uq is 
the height of the energy barrier. The crossover tempera- 
ture T* , at which both factors are of equal significance is 
defined by the condition B = Uo/kT, therefore 



2Nj 
kB(a ) 



1 -y-\/i--o 2 



(14) 



It follows from the eq. ( 14 ) that T* -> if a -> 



Rare-earth triangular clusters. The archetype 
of the noncollinear Ising model is a recently synthesized 
molecular dysprosium triangle 842|. The standard, based 
on the path integral, quasiclassical technique to handle a 



problem of macroscopic quantum tunneling 16 17 19 -21 
is only qualitatively suitable for rare-earth triangles due to 
low spin of the system. To get quantitative results we have 
to deal with the solutions of the Schroedinger equation in 
the relative Hilbert space. 

There are eight possible orderings of the ion spins in 
the rare-earth triangular cluster, which are Xi = I + ++}> 

X2 = | - + +), X3 = M h), X4 = I + + ->, and %n = 

Xn-4(+ <-> — ) for n = 5, 6, 7, 8, see eq. 

In the basis of the eight vectors 

■01=Xl) ^2 = (X2 + X3 + X4)/V3, 

"03 = (2X2 - X3 ~ X4)/\/6, 1p4 = (X3 - X4)/\/2, 

(V'n+4 are obtained as ipn by the replacement Xn ~> Xn+4 
for n = 1, 2, 3, 4), the matrix of the Hamiltonian in eq. ([2| 
in a zero magnetic field reads 



where 
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Fig. 2: The Rabi-type oscillations between states |r w ±3), the 
longer period corresponding to A = 0.003j spans over eight 
shorter periods corresponding to the doubled splitting A = 
O.OOBj, thus showing the cubical oscillation period dependence 
on the splitting. 



It is possible to obtain the eigenvalues and eigenstates 
of the Hamiltonian, = E n $ n (n — 1, 8). The rele- 
vant explicit expressions are given in table [l] Clearly, the 
eigenvectors $ n (n = 1, 8) are an orthonormal set. 

The eigenstate with the lowest energy is $2, which is 
the superposition of the states close to the states with 
|t = ±3), i.e. $2 = (ip+ + tp-) /v2, where we make use 
of ip± = \/2 (f + (-x) ■ ipx >5 + f-(-x) ■ i>2,i)- The expected 
values of the toroidal moment in these states are 



1 -x 
y/l-2x + 



4x 2 



(16) 



Supposed that initially the system is in the |V'+)-state, 
one can obtain the probability P(t) to find the system in 
the |^_)-state by the time moment of t 



8 / . 

^<$ n |^_)<V+|$„) exp --E n t 



(17) 



where \ip t ) = exp (-%Ht) 

The probability function P(t) oscillates (see fig.[2|, thus 
giving the clear evidence of the macroscopic quantum tun- 
neling of the toroidal moment in the system. If the level 
splitting A is small compared to the exchange constant j, 
the oscillation frequency depends on parameter x = A/j 
as v — ax 3 with a — 48 GHz, so for the typical value 
t = v- 1 - 1 ms 8| we have x ~ 0.003, i.e. A - 0.03 
cnr 1 if j = 10 cm . The expected value of the toroidal 
moment given by eq. ( j 1 6[ ) is very close to ±3, namely 
(f) w ±3(1 - x 2 /2) - ±(3 - 1 • 10- 5 ). 

Let us now consider the influence of magnetic and elec- 
tric fields on the MQT of the toroidal moment. The states 
of a rare-earth ion of the cluster in the crystal field are 



p-4 



Macroscopic quantum dynamics of toroidal moment in Ising-type rare-earth clusters 



Table 1: The eigenenergies e„ 
designations are used: e±(x) — - 



En/j and the eigenstates (n = 1, 



of the Hamiltonian in eq. (15 1, the following 
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i±f^M and/±(x) 



1 ± |£y, where d(a;) = Vl + 2a; + 4a; 2 and x = A/j. 
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|1) = \a), |2) = |6), and |fe), the energy levels are —A/2, 
A/2, and Wfc ~ 100 cm -1 respectively. In the general case, 
the projection of the Hamiltonian V = Hcf + Hz (with 
the Zceman part %z = SA*bHJ) onto the quasidoublet 
states gives Vij = Vy - Y,k VikV kj /W k , i,j = 1, 2. 

In the case of the magnetic field applied perpendicular 
to the plane of the ion triangle 



A 
2" 



qiH A 



A 



V 22 = J- Q2H Z 



V u = V 21 = 0, 
(18) 

where q ll2 = g2 f i%^ k ^(l,2\J y \k}(k\J v \l,2).^ 

The symmetry of the system is not broken in the per- 



pendicular field, which, as it follows from eq. (18), results 
in the renormalization of the splitting A only 

A -> A = A + ( qi - q 2 )H 2 . 

Moreover, in the case of Kramers (dysprosium) ions, 
such a field produces a small splitting A ~ g y H of the 
ground state due to the smallness of g y , thus removing the 
degeneration. Since the typical value of relaxation time 
is 1 ms M, the splitting A is estimated at 0.03 cm -1 . 
The factor g y then equals to 0.06 in a 1 T magnetic field. 
This way, the theory developed is valid not only for non- 
Kramers but also for Kramers Ising-type ions as well. 

An external static electric field perpendicular to the 
plane of the ion triangle has no influence on the consid- 
ered system, because the MEE properties of the rare-earth 
triangular clusters are caused by the field, directed in the 



plane of the ion triangle 10 



The current driven dynamics. — Let us consider 
now the interaction between the toroidal moment and an 
external electric current. To describe the interactio n w e 
put the term V — — jT into the Hamiltonian in eq. (15). 



The electric current j could be the displacement current 
J = 4tt W ■ If ^ e e l ec t r i c field strength E linearly depends 
on the time and is directed along the Z-axis, then V = 
- • To • t, where v = dEz/dt. The f-operator matrix has 
the diagonal form r = (+3, +1, +1, +1, —3, —1, —1, —1) in 
the basis of ip n (n — 1, 8). The Hamiltonian H + V of 



the system can also be diagonalized, the eigenvalues En 
and eigenvectors $^' ) (n = 1, 8) are not given here due 
to their unhandiness. The energies of the two low-lying 
levels could be approximated (for x « 1) in the vicinity 
of the avoided level crossing at j z = as 



E± (x, j z 



E T 




3T • j s 



(19) 



t(x) 



(20) 



where the crossing energy is Eq = —\j and the avoided 
level splitting is 5 = ^jx 3 . 

The toroidal moment is a conjugate variable against the 
electric current density and can be found at zero temper- 
ature as 

c dE_{x,j z ) 
4irT dj z 

The plot of the relative (equilibrium) dependence is given 
in fig. [3] (curve 1). It is seen, that the current changes the 
direction of the spin twisting to the opposite, thus revers- 
ing the toroidal moment of the system in a way of the 
relatively sharp jump of At = 6. The anapole moment 
reversal, which could be called a reanapolization, requires 
the currents of 10 7 A/cm 2 , which seems to be quite achiev- 
able in experiment. 

The tunneling processes when sweeping the current j z at 
a constant rate over an avoided energy level crossing can 
be treated in the frames of the Landau-Zener tunneling 
model. The probability P to change the state character- 
ized by the quantum number of the toroidal moment at 
the avoided level crossing is given by the expression 



P = 1 — exp 



t:6 2 



{A-K/c)h-T -\T l -T 2 \-\dj z /dt\ 



, (21) 



where t% = — 3 and t 2 = +3 are the toroidal moment 
quantum numbers of the avoided level crossing with the 
splitting S = §j£ 3 , and dj z /dt is the constant current 
sweeping rate. For the rate of dj z /dt ~ 10 13 A/(cm 2 -s) 
and x = 0.003 we have P 0.95. 

With the Landau-Zener-Stiickelberg model in mind, we 
can now start to understand qualitatively the hysteresis in 
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Current density, 10 A/cm 



Fig. 3: The zero temperature plot of toroidal (anapole) mo- 
ment projection onto the laboratory Z-axis vs the external cur- 
rent density for A = 0.003j: thicker curve (1) in equilibrium, 
the other curves in the case of the finite value of current sweep- 
ing rate 10 13 A/(cm 2 -s), (2) from r = +3 to r = -2.7 and (3) 
from r = —3 to r = +2.7. 



the system considered (see fig. [3]). Let us start at a suffi- 
ciently large negative current j z . At very low temperature, 
all molecules are in the |r = —3) ground state. When the 
current is come down to zero, all molecules will stay in the 
|r = —3) ground state. 

When passing the current over the avoided level crossing 
region at j z ps 0, there is Landau-Zener probability P (0 < 
P < 1) to tunnel from the \r = —3) to the |t = +3) state. 
The dependence of the toroidal moment on the current 
then reads as follows 

r(x,j z ) = P ■ dE +(*>M + (1 - P) • dE -(y*) (22) 



djz 



dj z 



shown in fig. [3] (curve 3). The toroidal moment undergoes 
the jump of At = 6P. So the average value of toroidal mo- 
ment t becomes less than 3, namely r = 2.7 for P = 0.95. 
The relevant mixed quantum state relaxes to the equilib- 
rium |t = +3) state due to interaction with environment 
(the 3 — > 1 process in fig. [3| . 

Sweeping the current in the opposite direction from the 
|r = +3) state, we arrive to the mixed state with toroidal 
moment average value of —2.7, see curve 2 in fig. [3] The 
state relaxes as above to the equilibrium |r = —3) state 
(the 2 — > 1 process in fig. |3| . This way, we will come to 
the closed hysteresis loop. 

Conclusion. — In the present work, it is shown for the 
first time that there exists a possibility of the macroscopic 
quantum tunneling of toroidal moment (spin chirality) in 
the molecular system based on Ising-type rare-earth ions. 
It is highly important that the ground state of such mag- 
netically ordered systems has zero magnetic moment and 
should be characterized by toroidal moment. This is the 
reason why the slow relaxation observed experimentally in 



such systems cannot be caused by the MQT of the mag- 
netic moment, but can possibly be accounted for the MQT 
of the toroidal moment instead. 
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